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Abstract 



Connected radio interferometers are sometimes used in the tied-array mode: 
signals from antenna elements are coherently added and the sum signal applied 
to a VLBI backend or pulsar processing machine. Usually there is no computer- 
controlled amplitude weighting in the existing radio interferometer facilities. 
Radio frequency interference (RFI) mitigation with phase-only adaptive beam- 



forming is proposed for this mode of observation. Small phase perturbations are 
introduced in each of the antenna's signal. The values of these perturbations 
are optimized in such a way that the signal from a radio source of interest is pre- 
served and RFI signals suppressed. An evolutionary programming algorithm 
is used for this task. Computer simulations, made for both one-dimensional 
and two-dimensional array set-ups, show considerable suppression of RFI and 
acceptable changes to the main array beam in the radio source direction. 



1 Introduction 

Suppression of radio frequency interference (RFI) with adaptive beamforming is 
widely used in radio astronomy, radar and telecommunications. The main idea be- 
hind many algorithms proposed for use in radio astronomy consists of weighting the 
outputs of array elements in such a way as to create zero values in the synthesized 
array pattern in the direction of RFI and to keep the signal of interest (SOI), the 
radio source to be observed, in the maximum of the main lobe without significant 
loss of gain [l~5l El [9]. During recent years there has been a growing interest in 
radio astronomy for applying these methods of RFI mitigation both to existing ra- 
dio telescopes and to future generation projects. There are several specific features 
of the large connected radio interferometers (RI) used in radio astronomy such as 
Westerbork Synthesis Radio Telescope (WSRT), Very Large Array (VLA) and Giant 
Metrewave Radio Telescope (GMRT) which make the straightforward application of 
this adaptive beam-forming different and difficult when compared to classic phased 
arrays: 

1. Connected RI are highly sparse arrays. 

2. Their main mode of operation is correlation processing. 

3. Direction of arrival (DOA) of a signal of interest is a known and time-dependent 
vector. 

4. There is no computer-controlled amplitude weighting in the existing RI backend 
hardware. 
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5. There is an auxiliary tied-array facility which is used during VLBI and pulsar 
observations. The mode of observation is similar to that of standard phase arrays: 
the signals from the antennas are added but without amplitude weighting, because 
the antennas of RI are identical. There is a phase-only control allowing coherent 
adding. The RI works as a "single dish". 

6. Noise-like radio source signals are usually much weaker than system noise 
(antenna + receiver) and RFI. 

Phase-only adaptive nulling is proposed for RFI mitigation during tied-array ob- 
servations. Small phase perturbations are introduced to the signals of every antenna. 
The values of these perturbations are optimized in such a way that the signal from 
SOI is preserved and the RFI signals suppressed. This techniques has been widely 
discussed [HI QUI EH El El 121 EI2] and is well suited to tied-array observations. 

2 Narrow-band model of SOI and RFI 

There are two approaches to adaptive beam-forming: narrow-band (complex weight- 
ing of amplitudes and phases) and wide-band (digital filtering, delay-tap weighting). 
The narrow-band approach will be used in the following text. 

Let us consider an equidistant M-element linear array. The M-dimensional array 
output vector X(9), as a function of an angle, i. e., the complex amplitude of the 
temporal signal x(t) = X(^)e j,27r ^ *, consists of the following components: 

N 

X{9) = S(9 ) + RFIn{O n ) + N sys (1) 

71=1 

where S(8q) is the signal vector corresponding to the plane wave coming from 
the direction 9o, RFI n (9 n ) is the nth RFI vector, coming from any direction 9 n , N sys 
is the system noise vector. These three components are uncorrelated. Vector S(9 ) 
depends on the incidence angle 9 Q of the plane wave, measured with respect to the 
normal to linear array 

S(9 ) = [l,e-* po ,...e- i ( M - r >' po ]' r (2) 

where phase shift <p = (27id/\) sin(0 o ), d is the spacing between array elements, A is 
the wavelength. The phase of the first antenna is chosen to be equal to 0. 

The beamformer, in general, consists of the complex weights w m e l ^ m , m = 1...M, 
which form the beamformer vector W 

W = [l,w 2 e^\...w M e iHl } T . (3) 

The output of the phased array is 

Y = W H X. (4) 

The beamformer should satisfy both following requirements: 

a) steering capability: the SOI is protected (W H S = g), for a prescribed direction 
9o, the response of the array is constant regardless of what values are assigned to the 
weights W; 
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b) the effects of RFI should be minimized. 

The minimum-variance distortionless response (MVDR) beamforming algorithm, 
subject to this constraint when g — 1, is proposed in order to minimize the variance 
of the beamformer output pQ . The solution for W in this case is 

Wmdvr = R~ X S (0 O ) [S (9 ) H R- l S (6 )r\ (5) 

where R is the correlation matrix of X. 

As was mentioned in the introduction, WSRT and other large radio astronomy 
arrays cannot use this algorithm in real time with the existing equipment because 
there are no amplitude control facilities. That is to say that "RFI nulling" is limited 
to phase- only control. 

3 Phase-only adaptive nulling 

Phase-only weights can be found to be the solution to the following system of nonlinear 
equations: 



M 



Real{J2 ei<t>mS m = M} (6) 

m=l 
M 

Imag{Y^e l<t>m S m = 0} (7) 

m=l 

M 

Real{J2 e i(/)m RFI m = 0} (8) 

m=l 
M 

Imag{J2 e i(t>m RFI m = 0}. (9) 



m=l 



e i<t>m are W eights in our phase-only case. The vector S is known and is de- 
termined by the SOI coordinates. The construction of the RFI vector requires the 
knowledge of RFI's DOA, which may be known beforehand or could be obtained 
from the observed correlation matrix i?y =< XiXj >,i,j = 1...M, because SOI is 
always much weaker than RFI. But it is necessary to have a special correlator for this 
purpose in order to follow all rapid scintillations of RFI which are usually averaged 
by the main radio interferometer correlator. So, in principle, the system of equations 
(6-9) can be solved and the phase corrections <p m introduced into the phase control 
system. The optimal solution (5) can be used as a zero approximation for m . The 
difficulties in implementing the solution of the system (6 - 9) are not mentioned here. 
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4 Total power detector output 

A more practical method of calculating the phase corrections (f) m in the tied array 
case is proposed here. The tied array total power detector output (TPDta) is 

TPD TA =< [ [x(t)] 2 dt>=TPD sig + TPD RFI + TPD N , (10) 
J o 

where the total power components TPD s i g ,TPD RF i,TPD]sf correspond to the 
SOI, RFI and system noise, respectively, < ... > means the statistical expectation. 
The mean value of the signal component is 



TPD sig (<j) m ) =< 



J Q |Xj cos [ 27r ^ ot + 27l (™- l)dsin(# )/A + m ]cftj >, (11) 



the mean value of the n-th RFI component is 



TPD RFIjn {(j) m ) =< A R Fi,n cos[27r f t + 2vr(m - l)d sm(9 RFI>n )/X + <p m }dt \ >, 

Jo U=i J 

(12) 

the mean value of the system noise component is constant. We assume also that 
TPD sig -C TPD R Fj^ n and the different RFI n are uncorrelated. 

Considering the TPD output as a function of M variables 4> m , the following cri- 
terium for a "good" vector $ = [<pi.... 4>m] T can be proposed: 



TPD stg ^) 
Y,Z=i TPD RFi,n{®) +TPD N 



C($) = ^ N Sl9 ^ max . (13) 



The denominator Yl,n=iTPD R Fi : n{$) + TPD N is the total TPD output under the 
assumption TPD sig « TPD RFI n . The numerator TPD sig (Q) can be calculated for 
each given $ and 9q (the DOA of SOI). Therefore, maximizing C(<I ) ) with a proper 
choice of $, a higher signal-to- RFI-plus- noise ratio at the tied array output can be 
achieved. 

This is a classic M-variable optimization problem which is difficult to solve by the 
common gradient methods because of the multimodality of C($) : there are many 
local (secondary) maximums and a searching algorithm will "get stuck" at one of 
them without finding the global maximum. 

Genetic algorithms (GA) search the solution for the set of variables through the use 
of simulated evolution, i.e., survival of the fittest strategy. In contrast to the gradient 
algorithms which are, in general, calculus-based algorithms, GA, first introduced by 
[8], exploits a guided random techniques during optimization procedure [3[TT1I1]. The 
multimodality problem is successfully overcome by this algorithm. 

A simplified block diagram of GA implementation in a radio interferometer is 
depicted in Figure 1. A phase control subsystem introduces a certain initial phase 
distribution $o corresponding to radio source coordinates and preliminary phase cal- 
ibration corrections. The output of the TPD is then continiously measured and used 
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to supply the GA program with the data (cost function samples) which monitor the 
performance of the tied array with respect to RFI. The GA uses these data to calcu- 
late new phases $ m with the aim to maximize C($). These new phases are introduced 
into the phase control subsystem after each iteration and a new value of the TPD 
output signal is used for the next step. Thus the feedback loop, phase control sub- 
system - TPD - GA, maintains the low value of TPD T a(&) and therefore the high 
value of the fitness function C(<&), i. e., the high signal-to- RFI ratio. 

5 Computer simulation 

Computer simulation was performed to illustrate the effectiveness of the phase-only 
nulling in RFI mitigation. 

First, a 14-element half- wavelength linear array was modelled. The SOI direction 
is equal to 0°, and there are two RFI signals: one at the angle —20.1° , and the 
other at the angle +10.015°. Figure 2 shows in logarithmic scale the quiescent (dash 
line) and adapted (solid line) array patterns. The significant suppression of RFI with 
the adapted pattern is clearly visible, while the quiescent pattern has the secondary 
lobe maximums at the RFI positions. Figure 3 shows the corresponding array phase 
distribution. 

Sparse 14-element array patterns are represented in Figures 4, 5, 6 and 7. The 
distance between the elements is 144m and the central frequency is 1420MHz. The 
lobes are much narrower than for the half-wavelength array, so the different figures 
are given to illustrate the result of phase-only nulling. Figure 4 (a linear scale pre- 
sentation of the pattern) is given here to illustrate the loss and distortions of the 
main lobe. This is more visible in the linear scale, whereas RFI suppression is better 
seen in the logarithmic scale. Figure. 5 - the quiescent and adapted patterns around 
angle 0°, logarithmic scale, Figures 6, 7 - the same patterns around the directions of 
RFI1 (DOA= +10.015°) and RFI2 (DOA= -20.1°). The corresponding array phase 
distribution is shown in Fig. 8. 

The subsequent figures illustrate this phase-only nulling for a two-dimensional 
planar array. 

Adaptive nulling was simulated for the half-wavelength array with 10x10 elements, 
the central frequency is equal to 1420MHz. Rectangular coordinates al and a2 are 
angles measured from the x and y axes, respectively, to the line from the array to the 
radio source; thus the SOI is at the zenith, with coordinates (90°, 90°). Coordinates 
of RFI were chosen so as to put them on the maximums of the secondary lobes, the 
values of RFI suppression are shown in the captions. The following sequence of figures 
is given: 

Figure 9: normalized (A(90°, 90°) = 1) quiescent array with indicated RFI posi- 
tions; 

Figure 10: array's pattern after adaptation; 

Figure 11: this array's phase distribution after adaptation. 



RFI mitigation 



6 



6 Conclusions 

1. Existing large radio interferometers (WSRT, VLA, GMRT) have only phase control 
facilities and the real-time adaptive nulling in the RFI direction should take this 
constraint into account . 

2. The total power detector at the tied- array output can be used for phase-only 
RFI mitigation as an indicator of the level of RFI. 

3. The Genetic Algorithm is a convenient tool for cost function maximization 
during the search for the optimal array phase distribution. 

4. Computer simulations show significant RFI mitigation for the sparse linear 
array in the narrow-band approximation (A///o << 1). 

5. Phase-only nulling can also be used for real-time RFI mitigation at the station's 
level in new projects such as ATA, LOFAR and SKA. 
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Figure 1: Block diagram of TPD adaptive phase control 




Figure 2: Quiescent (dash line) and adapted (solid line) 14-element half- wavelength 
array pattern, logarithmic scale; two nulls in the adapted pattern at —20.1° and 
+10.015°. 
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Figure 3: Phase perturbations corresponding to the adapted array pattern in Figure 
2. 




Figure 4: Quiescent (dash line) and adapted (solid line) 14-element array pattern, 
spacing=144m, central frequency=1420MHz, main beam, linear scale. 



RFI mitigation 



10 




1CT' ' ' 

-0.02 -0.015 -0.01 -0.005 0.005 0.01 0.015 0.02 

angle 



Figure 5: Quiescent (dash line) and adapted (solid line) 14-element array pattern, 
spacing=144m, central frequency=1420MHz, main beam, logarithmic scale. 




Figure 6: Quiescent (dash line) and adapted (solid line) 14-element array pattern, 
spacing=144m, central frequency=1420MHz, direction +10.015°, logarithmic scale. 
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Figure 7: Quiescent (dash line) and adapted (solid line) 14-element array pattern, 
spacing=144m, central frequency=1420MHz, direction —20.1°, logarithmic scale. 
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Figure 8: Phase perturbations corresponding to the adapted array pattern in Figures 
5, 6 and 7. 



RFI mitigation 



12 



before adaptation 




a2 



Figure 9: Two-dimensional lOxlO-element half-wavelength, quiescent array pattern, 
linear scale, central frequency= 1420MHz, RFI-1 at [45°, 90°], RFI-2 at [90°, 135°], 
directions of RFI coincide with the maximums of the sidelobes. 
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Figure 10: Two-dimensional lOxlO-element half-wavelength, adapted array pattern, 
central frequency= 1420MHz, RFI-1 at [45°, 90°], RFI-2 at [90°, 135°], linear scale; 
RFI-1 suppression=106.2dB, RFI-2 suppression=103.1dB. 
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Figure 11: Phase perturbations corresponding to the adapted array pattern in Fig. 
10: upper panel shows a 3D-presentation of the phase surface, lower panel shows the 
phase distribution in linear order. 



